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Summary. We have measured the two-dimensional kinematics and line-strength dis- 
tributions of 72 representative nearby early-type galaxies, out to approximately one 
effective radius, with our panoramic integral-field spectrograph SAURON. The resulting 
maps reveal a rich variety in kinematical structures and linestrength distributions, in- 
dicating that early-type galaxies are more complex systems than often assumed. We are 
building detailed dynamical models for these galaxies, to derive their intrinsic shape 
and dynamical structure, and to determine the mass of the supermassive central black 
hole. Here we focus on two examples, the compact elliptical M32 and the E3 galaxy 
NGC4365. These objects represent two extreme cases: M32 has very regular kinemat- 
ics which can be represented accurately by an axisymmetric model in which all stars 
rotate around the short a:xis, while NGC4365 is a triaxial galaxy with a prominent kine- 
matically decoupled core, with an inner core that rotates about an axis that is nearly 
perpendicular to the rotation axis of the main body of the galaxy. Our dynamical mod- 
els for these objects demonstrate that two-dimensional observations are essential for 
deriving the intrinsic orbital structure and dark matter content of galaxies. 



1 The SAURON project 

The formation and evolution of galaxies is one of the most fundamental research 
topics in astrophysics. A key question in this field is whether early-type galaxies 
form very early in the history of the universe or are gradually built up by mergers 
and the infall of smaller objects. The answer to this problem is closely tied 
to the distribution of intrinsic shapes, the internal dynamics and linestrength 
distributions, and the demography of supermassive central black holes. 

In the few past decades, it has become clear that ellipticals, lenticulars, and 
spiral bulges display a variety of velocity fields and linestrength distributions. 
Two-dimensional spectroscopy of stars and gas is essential when attempting to 
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derive information on the intrinsic structure. For this reason, we have built a 
panoramic integral-field spectrograph, SAURON (PP), which provides large-scale 
two-dimensional kinematic and linestrength maps in a single observation. 

We commissioned SAURON on the 4.2m William Herschel Telescope on La 
Palma in 1999. In low-resolution mode, the spectrograph combines a large field- 
of-view (33"x41") with a pixel size of 0'.'94. When the seeing conditions are good, 
the high-resolution mode, with a pixel size of 0'.'28, allows zooming in on galactic 
nuclei. SAURON observes in the spectral range of 4810-5340 A, which contains 
the gaseous emission lines H/? and [OIII] and [NI] , as well as a number of stellar 
absorption features (Mg&, Fe, H/3). The instrumental dispersion is ^100 km/s. 
Between 1999 and 2002, we have used SAURON to observe a carefully-selected 
representative sample of 72 ellipticals, lenticulars and Sa bulges, distributed 
over a range of magnitudes, ellipticities, morphologies and environments f|29j). 

We have finalized the data reduction, have accurately separated the emission- 
and absorption lines, have calibrated the line-strength measurements, and have 
in hand maps of the stellar and gaseous kinematics and linestrengths for all 48 
E and SO objects, with those for the spirals to follow soon. The maps reveal 
many examples of minor axis rotation, decoupled cores, central stellar disks, 
non-axisymmetric and counter-rotating gaseous disks, and unusual line-strength 
distributions ([21131). We have also developed new methods to spatially bin the 
data cubes to a given signal-to-noise and to quantify the maps with Fourier 
methods (|S1E|)- This allows accurate measurements of, e.g., the opening angle 
of the isovelocity contours and of the angle between the direction of the zero- 
velocity contour and the minor axis of the surface brightness distribution ([H]), 
enabling various statistical investigations of the entire sample of objects. 

2 Dynamical Models 

We are constructing detailed dynamical models which fit all kinematics and 
eventually even observations of the stellar line-strengths of the galaxies in the 
SAURON survey. We do this by means of Schwarzschild's (JHI) orbit superposition 
method, which was originally developed to reproduce theoretical density distri- 
butions (e.g., j^^^EJElEOI)) and was subsequently adapted to incorporate 
observed kinematic data in spherical and axisymmetric geometry f[Hll??l [milfip . 
We have implemented a number of further extensions including the ability to 
deal with a Multi-Gaussian Expansion of the surface brightness distribution 
(inilHlEI)- We have also shown that the large data sets that are provided by 
instruments such as SAURON can be modelled without any problems (|26|'). 

Recently, we completed the non-trivial extension to the software that al- 
lows inclusion of kinematic measurements in triaxial geometry f|27|). As in the 
axisymmetric case, observational effects such as pixel binning and point-spread- 
function convolution are taken into account. The chaotic orbits are dealt with 
in the 'standard' way (see 123), and the line-of-sight velocity profile is used to 
constrain the models. In the next two sections, we describe two applications in 
more detail, one in axisymmetry and the other for a triaxial intrinsic shape. 
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Fig. 1. Top panels: the SAURON kinematic maps for M32. From left-to-right: the mean 
velocity, velocity dispersion and Gauss-Hermite parameters /13 and /14, which measure 
the first and second order deviations of the line-of-sight velocity distribution from a 
Gaussian shape. Bottom panels: idem, but now for the best-fit axisymmetric dynamical 
model with /-band M/L = I.8M0/L0, Mbh = 2.5 x 10*^1/0, and i = 70°. 

3 Axisymmetric models for M32 

We applied our axisymmetric modeling software to the nearby compact E3 
galaxy M32 ( 26 ) . By complementing the SAURON maps (Figure ^ with high- 
resolution major axis stellar kinematics taken with STIS (^01)) the models are 
constrained at both small and large radii, which allows us to measure an accurate 
central black hole mass Mbh, stellar mass-to-light ratio M/L, and inclination i. 
The left panels of Figure |21 show the dependence of Ax^, which is a measure of 
the discrepancy between model and data, on Mbh, M/L (in solar units, for the 
J-band) and i. The inner three contours show the formal 1, 2 and Sir-confidence 
levels for a distribution with three degrees of freedom. The black hole mass and 
mass-to-light ratio arc constrained tightly at M, = 2.5 x lO^M© and M/L = 
I.SMq/ Lq, and the inclination is constrained to a value near 70° ±5°. The right 
panels of the same figure show similar contours, but now for a data-set consisting 
of the STIS-kinematics together with four slits extracted from the SAURDN-data. 
In this case the constraints on all three parameters, but most notably on the 
inclination, are much less stringent. This demonstrates that two-dimensional 
observations are essential to gain insight into the intrinsic structure of galaxies. 
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Fig. 2. Dynamical models for M32 f |26|'). The panels show contours of the goodness-of- 
fit parameter Ax^ as a function of the central black hole mass Mbh, the stellar mass- 
to-light ratio M/L and the inclination i. Each dot represents a specific axisymmetric 
dynamical model. The intrinsic fiattening q of the models is indicated in the lower-right 
corner of each panel. The models are constrained by STIS kinematics along the major 
axis f|lU|) together with two-dimensional observations obtained with SAURON in its 
high resolution mode ('|29|'). The inner three contours represent the formal 1, 2 and 3a- 
confidence levels for a distribution with three degrees of freedom. Left panels: model fits 
to a data set consisting of the STIS-data and the full SAURON field. Tight constraints are 
placed on the central black hole mass and mass-to-light ratio, as well as on the allowed 
range of inclinations. Right panels: the Ax^ for models that were constrained by four 
extracted slits from the 9" x 11" SAURON field (major and minor axis, and at ±45°, as 
in [lip and the STIS data. This shows that the traditional kinematic coverage provides 
almost no constraint on i, and that the resulting uncertainties on the inferred values 
of M/L and Mbh are correspondingly larger. 
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4 The triaxial galaxy NGC 4365 

The upper panels of Figure|21show the steUar kinematics in the central 30" x 60" 
of the giant elliptical galaxy NGC4365, derived from two SAURDN pointings ([7]). 
The velocity field clearly shows a prominent decoupled core in the inner 3" x 7" 
(cf. [23]). It has a rotation axis which lies 82° ± 2° away from that of the body of 
the galaxy, which rotates around its long axis. Such a structure is possible when 
the shape is intrinsically triaxial because of the presence of orbits that have net 
mean streaming around either the long or the short axis. 

The globular cluster system of NGC 4365 shows evidence for an intermediate 
age population (j^). The SAURON linestrength maps, however, indicate a pre- 
dominantly old stellar population suggesting that the observed kinematic 
structure may have been in place for over 12 Gyr and the galaxy is in stable 
triaxial equilibrium. We therefore applied our developed modeling software to 
this case, to investigate whether it is possible to reproduce all the kinematic data 
in detail, and to constrain M/ L and the intrinsic shape and orbital structure. 

We represented the observed surface brightness distribution of NGC 4365 by 
a Multi-Gaussian Expansion which accurately fits the observed radial variation 
of ellipticity, the boxyness of the isophotes, and the modest isophotal twisting. 
We derived the deprojected density by assuming that each of the constituent 
Gaussian components is stratified on similar concentric triaxial ellipsoids. The 
three Euler angles that specify the orientation of the ellipsoids can be chosen 
freely. For each choice, we computed a library of 4000 orbits, obtained from 20 
energy shells with 200 orbits each, covering the four major orbit families, and 
including orbits from minor families and chaotic orbits. As the spatial resolution 
of the SAURON measurements is modest, we did not consider the effect of a central 
black hole. The preliminary results indicate that the quality-of-fit parameter A'x^ 
varies quite significantly with M/L and the parameters defining the intrinsic 
shape. The lower panels of Figure O show the predictions of one model that fits 
the data well. This illustrates that the software works, and shows that NGC 
4365 is indeed consistent with a triaxial equilibrium shape. 

In principle, best-fit values of the shape parameters, the direction of observa- 
tion, and the mass-to-light ratio can be determined by a systematic investigation 
of the parameter space, just as was done for M32. For triaxial systems this is a 
very time-consuming effort, but a first-order guess of the galaxy parameters can 
be obtained by using other, simpler, schemes (see, e.g., [21I221I2H1)- Our detailed 
dynamical modeling software can then be used to explore this more restricted 
parameter range. Work along these lines is in progress, and will make it possible 
to deduce, e.g., the intrinsic properties of the kinematically decoupled cores seen 
in many of these systems. Inclusion of higher spatial resolution data will allow 
accurate measurement of the mass of the central black hole. 
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Fig. 3. Observations and dynamical models for the E3 galaxy NGC4365. Top panels: 
from left to right, the stellar velocity field, velocity dispersion, and Gauss-Hermite 
moments hs and hi, as observed with SAURON. The maps are based on two semi- 
overlapping pointings, sampled at O'.'S x O'.'S, and were constructed via a kinemetric 
expansion to provide the best representation of the data that is consistent with an 
intrinsically triaxial geometry (e.g., point-antisymmetry for the V and /13 maps, cf. 
0). The original maps can be found in (|7|). The amplitude of the velocity field is 
about 60 km/s, the peak velocity dispersion is 275 km/s, and the contours in the ha 
and /i4 maps range between ±0.10. The decoupled core measures 3" x 7". Bottom 
panels: idem, but now for a dynamical model with average intrinsic axis ratios p = 0.93 
and q = 0.69 (triaxiality parameter T — [1 — — g^) ~ 0.22), observed from a 

direction defined by the viewing angles = 85° and if — 15°. This model reproduces 
all the main characteristics of the observations ( 28 ). 
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5 Concluding remarks 

We have presented two examples of recent results from our program to con- 
struct detailed axisymmetric and triaxial dynamical models for galaxies in the 
SAURON representative survey of nearby ellipticals, lenticulars and Sa bulges. 
The panoramic SAURON observations tighten the constraints on the possible ori- 
entation of a galaxy considerably. The extension of the modeling software to 
triaxial shapes including kinematic constraints works, and that it will help us 
gain significant insight into the structure of early-type galaxies. 
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